#!/usr/bin/tcsh
# a script to check NAG routine E01AEF for a Chebyshev polynomial interpolant
# ifort e02a*.f e01ae*.f ~/lib/liblnag.a
# ifort e01aefe.f e01aeft.f e01aeut.f e01aevt.f e01aewt.f e01aext.f e01aezt.f e01aeyt.f e02afft.f e02ahzt.f e02akyt.f
# ifort e01aefe.f interpolants.f  e02afft.f e02ahzt.f e02akyt.f e02akzt.f p01abft.f p01abzt.f x01aaft.f x02ajft.f
# ifort e01aefe.f interpolants.f polyfits.f p01abft.f p01abzt.f x01aaft.f x02ajft.f x02amft.f x04aaft.f x04baft.f
# next line is for standard NAG e01aefe.f
# ifort e01aefe.f interpolants.f polyfits.f p0x0.f
# next line is for my test e01aefeMy.f
echo 'testing e01aef.f with e01aefeMy.f'
gfortran -ffixed-line-length-132 -g e01aefeMy.f interpolants.f polyfits.f p0x0.f
rm -f qq.r
./a.out < ./e01aefe.d
# ./a.out < ./e01aefe.d > & qq.r
# diff qq.r e01aefe.r
